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PREFACE 


This  paper  is  the  logical  continuation  of  the  research  started  with 
Report  ORA-63  -5 .  In  the  introduction  of  that  report  the  author  expressed 
his  hope  that  he  would  he  able  to  investigate  the  influence  of  such  phe¬ 
nomena  as  boundary  layers  and  shock  fronts  on  the  Doppler  shift.  The 
deeper  the  author  penetrated  into  those  problems  the  more  it  became  ob¬ 
vious  that  first  the  problem  of  the  influence  of  the  "duration  of  obser¬ 
vation"  on  the  Doppler  shift  needs  to  be  clarified.  This  is  the  time 
required  to  measure  the  frequency  shift.  The  present  report  is  an  analy¬ 
sis  of  the  influence  of  the  duration  of  observation  on  the  Doppler  shift. 
The  boundary  layer  problem  and  the  shock  problem  are  deferred. 

Frequent  references  are  made  to  Report  ORA-63 -5-  The  use  of  equa¬ 
tions  from  that  report  is  indicated  by  an  asterisk.  For  example,  (3*7)* 
means  equation  (3-7)  of  ORA-63  -5  • 
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ABSTRACT 


Formulae  are  developed  for  the  Doppler  shift  at  signals  that  either 
go  directly  from  the  emitter  E  to  the  missile  M  or  that  reach  the 
receiver  R  after  they  have  been  reflected  from  M.  E,  M,  and  R  may 
have  any  velocity  relative  to  the  heterogeneous  atmosphere. 

The  derivation  of  the  formulae  accounts  for  the  influence  of  the  time 
needed  for  the  frequency  measurement  on  the  Doppler  shift.  As  a  conse¬ 
quence,  the  refractive  index,  the  velocities  or  velocity  components  that 
enter  the  formulae,  are  mean  values  over  well  defined  intervals  in  space 
or  time. 

The  relation  between  the  error  in  the  velocity  measurement  and  the 
time  measurement  is  established. 

Keywords : 

Doppler  Effect 

Doppler  Radar 

Refractive  Index 


v 


ACKNOWLEDGMENT 


The  author  owes  a  debt  of  gratitude  to  Dr.  H. 
Mr.  H.  W.  Engelhardt  of  the  Computation  Division, 
Development  Center,  for  many  valuable  discussions 


R.  Lotze  and  to 
Air  Force  Missile 


vi 


TABLE  OF  CONTENTS 


Chapter  Page 

INTRODUCTION  .......  .  1 

'  I .  DESCRIPTION  OF  THE  PROBLEM .  3 

II.  THE  DURATION  OF  OBSERVATION . . .  9 

III.  THE  TWO  PEAKS  OF  VIBRATION  LIMITING  THE  DURATION  OF 

OBSERVATION  ........  .  13 

IV.  DOPPLER  SHIFT  FOR  STATIONARY  EMITTER .  17 

V.  MOVING  EMITTER  AND  RECEIVER .  27 


VI 


ACCURACY  CONSIDERATIONS 


31 


FINE  STRUCTURE  OF  THE  DOPPLER  SHIFT 


IN  HETEROGENEOUS  MEDIA 

INTRODUCTION 

The  method  employed  in  ORA-65-5  consisted  in  the  repeated  application 
of  the  principle  of  the  invariance  of  the  phase  function  under  the  Lorentz 
transformation  and  of  the  principle  of  the  constancy  of  the  frequency  in 
distinguished  frames  of  reference.  The  power,  usefulness,  and  mathemati¬ 
cal  elegance  of  that  method  need  no  further  demonstration.  Nevertheless, 
the  method  proves  rather  inadequate  for  the  treatment  of  the  following 
two  aspects  of  the  Doppler  shift: 

1.  The  influence  of  the  time  needed  to  measure  the  frequency  shift- 

2.  The  influence  of  such  disturbances  as  boundary  layers  or  shock 
fronts . 

To  account  for  those  genuinely  physical  details  the  method  of  ORA-63-5 
is  too  formal  and  too  mathematical.  What  is  needed  here  is  a  description 
of  the  space-time  geometry  of  the  light  and  missile  trajectories  during 
the  time  of  observation.  From  there  one  can  calculate  the  frequency  ratio 
in  terms  of  proper  mean  values  of  the  refractive  index  and  the  velocities . 


TABLE  I 


DEFINITIONS  OF  QUANTITIES  AND  SUBSCRIPTS 
E  =  emitter,  M  =  missile,  R  =  receiver,  T  =  atmosphere 
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CHAPTER  I 


DESCRIPTION  OF  THE  PROBLEM 


Throughout  this  report  we  shall  assume  that  the  atmosphere^  T  con¬ 
stitutes  an  inertial  body  of  reference.  By  this  we  mean: 

1 .  There  are  no  motions  between  parts  of  the  atmosphere .  With  this 
assumption  we  neglect  the  existence  of  boundary  layers  and  suppose  that 
the  light  trajectories  do  not  pass  through  the  shock  front  ahead  of  the 
missile  or  through  the  exhaust  plume  behind  the  missile . 

2.  The  atmosphere  is  at  rest  relative  to  the  earth. 

3 .  The  earth  itself  constitutes  an  inertial  body  of  reference . 

Of  course.,  the  last  assumption  is  not  strictly  realistic.  The  in¬ 
fluence  of  the  non-linear  motion  of  the  earth  on  the  frequency  shift  is 
probably  always  ne glee table,  but  the  influence  of  the  gravitational  field 
might  become  perceivable  under  certain  conditions.  This  problem  must  be 
deferred  to  later  investigations.  However,  we  can  safely  say  that  the 
gravitational  frequency  shift  destroys  itself  if  emitter  and  receiver  are 
on  the  same  gravitational  potential.  This  is,  for  instance,  the  case 
when  the  signal  is  emitted  on  ground  level,  reflected  from  the  missile, 
and  received  on  ground  level. 

The  used  electromagnetic  radiation  is  always  assumed  to  be  as  mono¬ 
chromatic  as  possible.  Since  we  consider  only  very  narrow  sections  of 

^We  use  the  symbol  T  for  the  atmosphere  in  order  to  be  consistent  with 
ORA-63 -5  where  it  means  "transparent  medium." 


the  waves  (called  beams),  the  waves  may  be  considered  as  plane.  A  mono¬ 
chromatic,  plane  wave  can  be  described  by  a  wave  number  "vector"  1c  and 
an.'  angular  frequency  to.  Like  the  ordinary  (classical)  velocity  v,  Tc 
is  not  a  vector  under  the  In  rent z  transformation  (LT).  And  like  the  re¬ 
fractive  index  n,  the  frequency  m  is  not  a  scalar  under  the  LT.  All 
these  quantities  depend  on  the  frame  of  reference.  Further,  k,  n,  and 
(»)  depend  on  the  location  in  the  atmosphere.  Finally,  the  velocities  must 
be  specified  as  to  the  objects  (emitters,  missiles,  mirrors,  receivers) 
that  have  these  velocities.  Hence,  all  of  the  four  quantities  k,  to,  n, 
and  v  need  two  specifications;  one  for  the  frame  of  reference  and  one 
for  the  location  or  the  object  that  has  the  respective  velocity. 

However,  in  this  report  (as  in  ORA- 63 -5 )  the  velocities  and  the  re¬ 
fractive  index  refer  always  to  the  atmosphere  T  so  that  the  subscript 
for  the  frame  of  reference  becomes  dispensable  for  velocities  and  refrac¬ 
tive  index.  We  therefore  adopt  here  the  same  subscript  notation  as  in 
ORA-63 -5 .  This  is  outlined  in  more  detail  in  Table  I,  page  2.  There, 
we  introduce  at  the  same  time  some  other  quantities . 

Let  now  kte,  ktm,  ktm,  ktr,  ve,  vm,  vr  denote  the  absolute 
amounts  of  the  respective  "vectors"  as  defined  in  Table  I.  We  then 
define  the  angles  te,  tm,  f  as  observed  from  T  at  the  re¬ 

spective  locations,  and  the  "normal  components"  ven,  vmn,  vmri,  Tm 
of  the  respective  velocities  as  follows; 


v  =  k,  v  cos  t  =  k.  v 
e  te  e  e  te  en 


(1.1) 
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(1.2) 


vm  "  ktm  vm  cos 


vm  23  ktm  vm  cos  fm 


v  =  k.  v  cos  t 
r  tr  r  r 


ktm  vmn 

=  ktm  ^nn 

=  k,  v 
tr  m 


(1.3) 

(l.fc) 


Figure  1  is  a  perspective  illustration  of  these  definitions.  The  word 
"perspective"  may  remind  the  reader  that  the  "vectors"  do  not  all  lie  in 
one  plane . 

The  points  E,  M,  and  R  are  the  instantaneous  (not  simultaneous'.) 
positions  of  emitter,  missile,  and  receiver.  The  directed  curves  S  and 
S  are  the  light  trajectories  from  E  to  M  and  from  M  to  R  relative 
to  T. 

We  shall  always  assume  that  the  refractive  index  -  as  far  as  its  geo¬ 
metrical  dependency  is  concerned  -  depends  only  on  the  distance  from  the 
center  of  the  earth  C  and  not  on  the  azimuth  angle.  An  immediate  conse¬ 
quence  of  this  assumption  is  that  all  light  trajectories  are  plane  curves . 
A  light  trajectory  going  from  E  to  M,  for  instance,  lies  in  the  plane 
EMC. 

For  our  later  purposes  we  define  now  four  cases  of  Doppler  shift  and 
list  the  formulae  as  they  have  "been  derived  in  ORA-63 -5  » 

CASE  Is  Emitter  E  stationary  on  the  ground  (ve  =  0)  and  missile  M 
moving  with  vm.  Signal  emitted  hy  E  and  received  by  M.  The 
corresponding  formula  is  obtained  from  the  more  general  formula 
(3.7)*  if  one  specializes  on  ve  =  0  and  replaces  the  subscript 
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r  by  m.  This  yields 


with 


_ 


%  P] 


mn 


w. 


ee 


V* 


(1.5) 


-in 


m 


„  vra 


(1.6) 


-inn 


-  Pm  cos  t, 


m 


(1.7) 


->  * 

Case  II:  Again  ve  =  0  and  vm  0,  but  signal  emitted  by  E,  reflected 

from  M  and  received  by  E.  This  is  the  case  for  which  formula 
(4.8)*  has  been  developed.  It  is 


wee.1+niP: 


w. 


ee 


mn 
^  ^mn 


(1.8) 


Case  III:  vg  0^  vm  4  signal  emitted  by  E  and  received  by  M. 
Slightly  altering  formula  (5.7)*  yields 
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with 


(1.10) 


"en 


=  Isa  =  p 

c 


e  cos  ^e 


(1.11) 


Case  IV:  This  is  the  most  general  case,  with  ^  0,  v„tO,  vr  =f  0, 

where  the  signal  is  emitted  by  E,  reflected  from  M  and  re¬ 
ceived  by  R.  For  this  ease  refer  to  formula  (4.6)*.  It  is: 


urr 


(l  ~  ^  ^mn^  (l  ~  ne  Pen) 
(l  -  nm  Pmn)  (l  “  nr  3rn) 


(1.12) 


with 


(1.13) 


Pm  =  Isa  =  pr  cos  tT  (i.i4) 

c 

It  is  now  our  aim  to  rederive  these  formulae  in  such  a  way  that  we 
can  account  for  the  influence  of  the  "time  of  observation."  This  will  at 
the  same  time  provide  for  a  more  detailed  description  of  such  quantities 
as  ne,  ve,  ven,  te,  tm,  which,  as  of  now,  are  "local"  or  "in¬ 

stantaneous  values."  All  these  quantities  will  be  described  as  mean  val¬ 
ues  over  well-defined  intervals  in  space  or  time. 
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CHAPTER  II 


THE  DURATION  OF  OBSERVATION 


Every  frequency  measurement  needs  a  minimum  of  time If  the  frequency 
is  supposed  to  be  constant  in  the  time,,  the  minimum  time  required  to  meas¬ 
ure  it  is  the  period  of  vibration.  However*  to  enhance  the  accuracy  with 
which  the  frequency  is  to  be  determined  one  better  observes  over  a  much 
longer  time  interval  At  and  counts  the  number  of  vibrations.  This  time 
interval  At  which*  within  certain  limits  can  be  selected  at  liberty,  is 
the  "duration  of  observation." 

It  is  now  our  aim  to  determine  the  duration  of  observation  which  is 
necessary  to  measure  the  Doppler  shift.  It  is  sufficient  if  we  confine 
our  consideration  in  this  chapter  to  the  technologically  most  important 
case,  which  is  Case  II  as  described  in  Chapter  I. 

The  techniques  of  measurement  are  always  more  or  less  interference 
methods .  A  part  of  the  original  wave  is  branched  off  before  emission  and 
"stored"  in  a  time  delay  device .  The  reflected  wave  is  grossly  amplified 
and  then  superimposed,  on  the  stored  wave.  If  A*  A*  <fb  and  <j>  denote 
the  amplitudes  and  the  phase  functions  of  the  original  and  the  reflected 
wave,  the  superposition  can  be  mathematically  described  by 


A  sin  <j>  +  A  sin  <i>  « 


A  +  A 


$  +  j>  ^ 


similar  statement  does  not  hold*  for  instance*  for  length  or  mass  meas¬ 
urement  . 
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If  one  then  filters  out  the  term  with  the  amplitude  (A  -  A)  or  if 
one  amplifies  such  that  A  -  A  =  0,  the  superposition  results  in  one 
single  modulated  wave  with  the  angular  modulation  frequency 


(2.1) 


In  order  to  measure  ui*  one  has  to  observe  at  least  one  period  of 
modulation  which  is 


T*  .  £2  =  - - Zl -  (2.2) 

l“ee  "  “eel  - 

We  call  t*  the  "minimum  duration  of  observation."  In  order  to  re¬ 
duce  the  error  of  measurement  one  may  observe  during  a  much  longer  period 
of  time  At.  However,  if  then  w*  and  T*  are  not  constant,  one  meas¬ 
ures  their  mean  value  over  the  period  At.  If  we  denote  these  mean  values 
by  (w*)  and  (t*)  and  if 

oil  (2.3) 

is  any  real  number  not  smaller  than  unity,  we  may  write 

At  =  a  (t*)  =  -  «  — — -  (2.4) 

<«*)  <|«tee-«ee|) 

This  result  shows  that  there  are  two  antagonistic  interests.  The 
first  one  calls  for  a  high  accuracy  in  measuring  («*) .  To  achieve  this 
one  would  have  to  make  At  and  a  great.  The  second  interest  is  to 


10 


obtain  information  about  (0*  which  is  as  detailed  as  possible.  This,  in 
turn,  calls  for  making  At  and  a  small.  Consequently,  one  has  to  seek 
a  compromise  between  these  two  interests;  that  is,  one  has  to  select  the 
number  a  or  the  duration  of  observation  At  within  certain  ranges.  In 
the  last  chapter  of  this  report  we  shall  return  to  this  interesting  subject. 

It  is  now  our  concern  to  estimate  the  order  of  magnitude  of  the  minimum 
time  of  observation  T*.  To  this  end  we  employ  formula  (1.8)  which  has  been 
derived  in  ORA-.63 -5 »  The  fact  that  we  disregarded  there  the  duration  of 
observation  as  well  as  the  boundary  layer  certainly  does  not  affect  the  mag¬ 
nitude  relations.  Slightly  rearranging  equation  (1.8)  yields 

<«>ee  "  “ee  “  Pjnn  (“ee  +  wee)  (2-5) 

For  all  technically  achievable  missile  velocities,  wee  and  u>ee  are 
very  close.  Further,  nm  is  very  close  to  unity.  Hence  we  obtain  from 
equations  (2.1)  and  (2.5) 

=  uee  I  Pmni  =  wee  Pm  I cos  $ml  (2.6) 

However,  we  need  (w*)  rather  than  w*.  Therefore,  we  formally  introduce 
the  mean  value  (pmn )  of  p^  but  leave  its  exact  definition  to  Chapter 
IV.  Hence,  we  replace  equation  (2.6)  by 

<«*)  =  wee  !<Pmn>l  s  wee  I K  cos  tm> |  (2.7) 

Substituting  this  into  equation  (2,k)  yields  then 

At  -  0  (t*>  -  =  JUSS-  (2.8) 

wee  I  ^mn^  I  I  ^Pmn^  I 
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Here,  Tee  is  the  period  of  vibration  of  the  original  signal  as  observed 
by  and  at  E  (see  also  Table  I). 

Equation  (2.8)  shows  that  the  duration  of  observation  is  far  greater 
than  the  period  of  vibration.  To  have  an  intuitive  example  let  us  assume 
a  typical  radar  frequency  of  some  109  cps  which  corresponds  to  a  period 
of  vibration 


Tee  ~  -^-0  s 

If  we  then  assume 


vm  -  3000  | 

and  exclude  extreme  cases  with  |  v^l  «  vm,  we  have  also 

I vmnl  «  3000  | 

With  c  =  ).0  X  10®  ~  this  yields 

IfW  - 10'5 

It  then  follows  from  equation  (2.8) 

<T*>  «  10"4  S 

which  means 

At  1  T*  »  Tee 
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CHAPTER  III 


THE  TWO  PEAKS  OF  VIBRATION 
LIMITING  THE  DURATION  OF  OBSERVATION 


As  we  mentioned  in  the  introduction,  we  shall  derive  the  Doppler  shift 
formulae  for  various  cases  in  a  certain  geometrical  manner.  This  method 
is  rather  laborious,  so  that  we  limit  its  detailed  outline  to  Case  I  of 
Chapter  I.  There,  E  is  at  rest  relative  to  T.  Hence,  clocks  that  rest 
in  E  (E-clocks)  and  clocks  that  rest  in  T  (T-clocks)  measure  the  same 
times .  If  we  do  not  specify  otherwise  the  times  in  this  and  the  next  chap¬ 
ter  are  measured  by  E-clocks  or  T-clocks . 

We  then  consider  a  certain  peak  of  vibration  of  the  monochromatic  radi¬ 
ation  emitted  by  E  and  received  by  M.  We  call  this  peak  the  "zero-peak." 
It  may  leave  E  at  the  time  teo  and  arrive  at  M  at  the  time  t^.  Its 
travel  time  is  then 


At 


mo 


eo 


(3-1) 


We  then  assume  that  the  number  a  in  equation  (2.8)  is  selected  such 
that  the  number 


p  s  At_  ^  (5.2) 

Tee  | (Pjm) | 

is  an  integer.  This  assumption  serves  only  to  make  the  derivations  of 
this  and  the  next  chapter  as  transparent  as  possible  and  is  in  no  way 
necessary  nor  even  essential.  In  particular,  this  assumption  has  no 
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impact  on  the  technique  of  measurement  even  though  a  is  generally  no 
longer  an  integer  if  P  is  an  integer  since  — -  is  generally  not 

KfWI 

an  integer.  When  performing  actual  measurements  one  is  not  bound  by  the 
assumptions  someone  made  in  his  theoretical  derivations. 

In  all  practical  cases,  P  is  a  very  great  number.  Under  our  spe¬ 
cial  assumption  it  is  now  a  very  great  integer.  Hence,  At  is  a  very 
great  multiple  of  Tee .  A  consequence  of  this  is  that  at  the  time 

tep  =  teo  +  At  (3.3) 

another  peak  of  vibration  leaves  E.  We  call  it  the  "P-peak."  The  P-peak 
is  not  consecutive  to  the  zero-peak  but  separated  from  it  by  a  very  great 
number  (P-l)  of  peaks.  Let  the  P-peak  arrive  at  M  at  the  time  t^. 
Its  travel  time  is  then 

Atp  “  ^mP  "  ^eP  (3*^) 

The  period  of  vibration  as  observed  by  and  at  E  is  simply 


T 


ee 


(3.5) 


However,  the  period  of  vibration  as  observed  by  and  at  M  is  not 
presented  by  the  analogous  expression,  because  it  is  measured  by  M-clocks 
while  tjjjQ  and  t^  are  measured  by  E-clocks .  Let  dTm  denote  the 
element  of  the  eigen  time^  for  bodies  that  are  at  rest  relative  to  M. 

’"’Here  only  one  subscript  makes  sense .  It  is  the  one  that  indicates  the 
frame  of  reference. 


l4 


If  then  dt  denotes  the  time  elements  of  the  E-clocks  it  holds  according 
to  the  special  theory  of  relativity 


drm  -  dt  yTT  f£(tj 

We  wrote  Pm(t)  in  order  to  indicate  that  it  may  vary  with  the  time. 

The  period  of  vibration  as  observed  by  and  at  M  becomes  then 

p  "^mP  .  ■■■  - - -  -y 

rmm  =  f  J  V 1  "  P®(t)  dt  (5'6) 

mo 

By  using  the  mean  value  theorem  we  may  write 

'am”  (-/win  ist/V  (5.7) 

with 


<VX "  ^ 


dt 


(3.8) 


Tee  and  Tjjjjjj  are  both  eigen  times  in  Minkowski’s  sense.  Therefore,  their 
ratio  is  the  inverse  of  the  frequency  ratio.  Hence,  we  have 


Tmm  Mee  „  ^  /l  ~  p2~'  )  ^  ~  ‘tm° 
Tee  '•^mm  »  m  ^eP  “  ^eo 


(3-9) 
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If  we  now  combine  equations  (3.1),  (3>*0>  (3  •  5 )  >  (3  •  T )  ^  and  (3  *9)  we  ob¬ 
tain 


ML 


mm 


(s5. 


ee 


<aA^> 


{ 


At, 


P  ~  ^o  1 


^mo  ^ 


(3.10) 


This  is  the 
the  ratio 


basic  equation  we  shall  use  in  the  next  chapter  to  calculate 
0Jmm 


'/ 


w  as  a  function  of  the  missile  velocity, 
ee 
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CHAPTER  IV 


DOPPLER  SHIFT  FOR  STATIONARY  EMITTER 


Equation  (3  .10)  at  the  end  of  the  last  chapter  sets  the  task  for  this 
chapter;  to  express  the  difference  Atp  -  At0  of  the  two  travel  times  in 
terms  of  the  missile  velocity.  This  is  a  geometrical  task  which  is  greatly 
facilitated  by  using  'T  as  frame  of  reference.  We  hope  that  this  will  be¬ 
come  clear  in  the  course  of  this  chapter;  however,  we  would  like  to  point 
out  now  that  employing  T  as  body  of  reference  is  not  a  necessity,  but 
rather  a  convenience.  On  the  other  hand,  a  convenience  can  amount  to  a 
necessity  when  the  difficulties  that  arise  from  disregarding  it  become  un- 
surmountable .  This  best  describes  the  present  situation. 

In  the  preceding  chapter  we  introduced  the  zero-peak  and  the  P-peak. 

In  this  chapter  we  have  to  consider  their  trajectories.  Now,  the  tra¬ 
jectory  of  one  single  peak  of  vibration  might  be  a  questionable  concept 
since  electrodynamics  as  well  as  classical  mechanics  or  thermodynamics  are 
maeroseopical  theories  and  as  such  statistical  in  nature.  Hence,  the  con¬ 
cept  "trajectory"  might  be  legitimate  only  if  it  refers  to  a  sufficiently 
great  number  of  waves .  However,  this  does  not  constitute  a  serious  prob¬ 
lem  in  the  present  case  since  the  number  P  of  vibrations  during  the  time 
of  observation  is  very  great  as  compared  with  unity,  (in  Chapter  VI  we 
shall  show  that  it  is  of  the  orders  107,  10s,  or  109.)  Consequently,  it 
is  always  possible  to  select  two  further  numbers  Px  and  P2  such  that 


1  «  Px  «  P 
1  «  P2  «  P 


(*.D 
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Figure  2.  Missile  trajectory,  light  trajectories,  and  phase  Surfaces 
with  respect  to  T  as  body  of  reference .  (Perspective 
picture) 
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Instead  of  considering  the  zero-peak  and  P-peak  one  may  then  consider 
the  first  and  the  last  P2  waves .  For  these  two  wave-groups  the  con¬ 

cept  of  a  trajectory  certainly  makes  sense.  On  the  other  hand,  the  number 
(P  -  (P]_  +  P2)}  of  waves  or  peaks  that  separate  the  two  wave-groups  is 
still  approximately  equal  to  P.  Yet  this  is  the  only  essential  we  shall 
use  as  far  as  the  zero-peak  and  the  P-peak  are  concerned.  Figure  2  is  a 
grossly  distorted,  perspective  picture  of  the  missile  trajectory,  the  tra¬ 
jectories  SQ  and  Sp  of  the  zero-peak  and  the  P-peak,  and  of  two  auxil¬ 
iary  light  trajectories  Sx  and  S2.  The  function  of  these  two  trajectories 
will  be  seen  later  in  this  chapter. 

Of  course,  the  light  trajectories  and  the  missile  trajectory  do  not 
generally  lie  in  one  plane.  Generally,  the  missile  trajectory  does  not 
lie  in  a  plane  at  all. 

The  points  Mq,  Mx,  M2,  Mp  denote  the  positions  of  the  missile  at 
the  times  tmo,  tmi,  t^,  tmp  (all  measured  by  E-clocks  or  T-clocks), 

The  curves  4>o,  4)1,  4>2  symbolize  surfaces  of  constant  phase  with  respect 
to  E  as  origin  of  radiation.  The  intersections  of  the  surfaces  4*0,  4^, 

4>2  with  the  curve  Sp  are  the  points  Aop,  Aip,  Asp.  The  intersection 
of  4>i  with  S2  is  the  point  A12. 

In  Figure  2  all  geometrical  quantities  like  trajectories,  surfaces, 
points,  and  angles  refer  to  T.  This  results  in  two  advantages  which  we 
shall  now  discuss  by  assuming  temporarily  that  T  as  frame  of  reference 
is  replaced  by  Z*  relative  to  which  T  is  in  motion.  We  indicate  this 
by  adding  an  asterisk  to  those  geometrical  quantities  we  deal  with. 
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The  heterogeneity  of  T  has  then  the  consequence  that  the. refractive 
index  n*  is  an  explicit  function  of  the  time.  If  we  then  consider  the 
two  light  trajectories  from  E*  to  along  SQ*  and  from  E*  to 

Aqp*  along  Sp*  we  can  make  the  .following  statement: 

Identical  waves  or  peaks  need  equal  times  for  both 
trajectories  because  and  AQp*  lie  on  one  sur¬ 

face  of  constant  phase,  but  different  waves  or  peaks 
need  different  times  even  for  one  and  the  same  tra¬ 
jectory  because  the  refractive  index  has  changed  in 
the  meantime . 

Only  if  T  itself  serves  as  frame  of  reference  need  different  peaks  or 
waves  take  the  same  time  for  the  same  trajectory.  If  we  now  return  to 
T  as  frame  of  reference,  At0  is  the  time  the  zero-peak  need  from  E 
to  Mq  as  well  as  from  E  to  A^.  On  the  other  hand,  the  P-peak  needs, 
from  E  to  Aop,  the  same  time  as  the  zero -peak.  If  we  denote  the  time 
of  the  P-peak  from  E  to  Agp  by  AtQp,  we  have 

At0  =  Atop  (4.2) 

This  is  the  first  advantage  of  employing  T  as  frame  of  reference. 

The  second  advantage  is  that  n  does  not  depend  on  the  direction, 
while  n*  does  because  of  the  drag  effect.  Hence,  T  is  isotropical 
only  with  respect  to  itself  but  not  with  respect  to  E*.  However,  in 
such  frames  of  reference  the  light  trajectories  and  the  surfaces  of  con¬ 
stant  phase  are  not  orthogonal.  Hence,  the  orthogonality  of  the  light 
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trajectories  and  the  phase  surfaces  in  Figure  2  is  the  second  advantage 
of  T  as  frame  of  reference. 

These  two  advantages  will  greatly  facilitate  the  following  deriva¬ 
tions.  Moreover,  it  is  necessary  for  the  understanding  of  the  entire 
complex  of  problems  that  one  is  conscious  of  these  facts. 

We  have  now  to  calculate  the  difference  Atp  -  At0  of  the  travel 
times  occurring  in  the  basic  equation  (3  .10).  The  travel  times  can  be 
expressed  by  line  integrals  along  the  respective  light  trajectories .  If 
dSQ,  <3%,  dS2,  and  dSp  denote  the  line  element  on  S0,  S1#  S2, 

and  Sp  we  may  write 


At0  = 


1 

c 


(*.3) 


(4.4) 


(*.5) 


M2 

At2  s  i  J  n2  dS2  (4.6) 

E 
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n2  denote  the  refractive  index  as  function  of  position 


nQ  f  np  f  nx y 

along  Sq  y  Sp  y  y  Sg  o 

By  construction  we  have  further 


1 

c 


(4-7) 


1 

c 


M2 

a12 


dSp 


If  we  now  Build  the  difference 


(4.8) 


Atp  -  At0 

and  observe  equations  (4.3),  (4.4),  (4.7)  and  Figure  2  we  find 


Atp 


dSp 


(4.9) 


Thus  we  have  to  express  this  line  integral  in  terms  of  the  missile  veloc¬ 
ity.  To  facilitate  this  we  introduced  the  auxiliary  light  trajectories 
Sx  and  S2 .  We  assume  that  the  time  difference 


dt  =  t2  -  tx  (4.10) 

is  sufficiently  small  so  that  the  triangle  MXA12M2  can  be  considered  as 
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composed  of  straight  lines.  Since  the  angle  at  A12  is  90  degrees  by- 
construction  we  have  then 


A12M2  =  MxM2  cos  fm  (4.11) 

where  tm  is  the  angle  between  the  two  "vectors"  and  vm  at  the 

position  under  consideration  and  as  introduced  in  Chapter  I. 

On  the  other  hand,  the  distance  MXM2  is  the  distance  the  missile 
has  traveled  in  the  time  dt  -  t2  -  tx .  Consequently,  we  have 

MiMg  =  vm(t)  dt 

If  we  substitute  this  into  equation  (4.1l)  we  obtain 

-  vm(t)  cos  tm(t)  dt  (4.12) 

That  (t2  -  tx)  is  infinitesimal  makes  the  integral  signs  in  equation 
(4.8)  dispensable.  Moreover,  the  line  element  dS2  can  be  replaced  by 
the  distance  A12M2.  If  we  then  write  n^  instead  of  n2  expressing 
that  this  is  the  value  of  n  in  the  immediate  neighborhood  of  M,  equa¬ 
tion  (4.8)  becomes 


%(■&)  vm(t)  cos  tm(t)  ==  np  dSp 

Integrating  the  left-hand  side  of  this  equation  from  t^  to  t^  corre¬ 
sponds  to  integrating  the  right-hand  side  from  Aq p  to  Mp.  This  yields 


r  ^mP  P 

~  /  nm(t)  vm(t)  cos  tm(t)  dt  =  i  /  np 


asr 


“mo 


Aop 
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If  we  then  substitute  this  into  equation  (4.9)  we  obtain 


tmP 


At0  -  Atn  = 


J  nm(t)  vm(t)  cos  ♦m(t)  dt 


^mo 


This  is  to  be  substituted  into  the  basic  equation  (3.10)  of  the  preceding 
chapter.  The  result  is 


mm  _  1 


.w 


Uee  (^/i  -  &)  L  (V  -  W 


J  nm(t)  Pm(t)  cos  tm(t)  dt  "j- 

(4.13) 


mo 


It  now  suggests  itself  to  apply  again  the  mean  value  theorem  and  to  write 


^mP 


Pm  cos  *m)  =  f  ^(t)  vm(t)  cos  tm(t)  dt  (4.l4) 

^  W 


^mP  "  tgn 


which  leads  to 


.  1  -  <°m  Pm  cos  »m> 

Wee  <V1  “ 

This  result  replaces  the  former  result  (l»5)»  It  has  the  same  alge¬ 
braic  structure  but  the  advantage  that  the  occurring  combinations  of  the 

l 

refractive  index  and  the  missile  velocity  are  now  mean  values  over  well- 
defined  time  or  space  intervals . 
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We  defer  all  further  discussions  of  this  result  to  Chapter  VI  and 
proceed  immediately  to  Case  II  of  Chapter  I  where  the  signal  is  reflected 
from  M  and  then  received  by  E.  This  complication  does  not  change  the 
geometry,  so  that  we  need  not  go  into  any  details  but  can  confine  ourselves 
to  listing  the  final  result: 

“ee  _  1  +  <%  cos  »m>  (4.l6) 

5  1  -  (n  B  cos  t  ) 

ee  'mm  m' 

This  equation  replaces  the  former  result  (1.8). 
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Figure  5-  Moving  emitter  and  moving  missile. 
(Perspective  picture) 
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CHAPTER  V 


MOVING  EMITTER  AND  RECEIVER 

We  turn  now  to  Cases  III  and  IV  of  Chapter  I.  We  need  not  go  through 
the  whole  analysis  but  may  confine  ourselves  to  those  special  features 
that  require  some  additional  thoughts .  This  concerns  primarily  the  geom¬ 
etry  which  is  pictured  in  Figure  3  for  Case  III. 

The  body  of  reference  is  T.  The  times  teo,  tep,  tmo,  t^,  tmo, 

t  t  ,  t  „  are  measured  by  T-clocks.  These  are  the  times  when  the 
mP;  ro'  rP  J 

zero-peak,  the  P-peak,  the  reflected  zero-peak,  or  the  reflected  P-peak 
leave  E,  arrive  at  M,  leave  M,  and  arrive  at  R,  respectively.  To 
these  events  correspond  the  positions  EQ,  Ep,  Mp,  RQ,  Rp  of 

emitter,  missile,  and  receiver.  (Rq  and  Rp  are  not  shown  in  Figure  3-) 
In  Figure  3>  SQ  and  Sp  are  the  light  trajectories  of  the  two  peaks 
from  E  to  M.  S  is  an  auxiliary  light  trajectory.  4>0  is  a  surface 
of  constant  phase  with  Ep  as  origin,  'j’p  is  a  surface  of  constant  phase 
with  Mq  as  imagined  origin  of  radiation.  It  then  holds  by  definition 


AdP  ^3 
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From  these  two  equations  follows 


dSp 


E 


P 


The  travel  times  for  the  two  peaks  are 


* 


(5-1) 

* 


(5-2) 


Hd  Mp 

Atp  =  —  J  np  <3Sp  »  —  J'  np  dSp  +  i  J  np  dSp  (5-3) 

Ep  Ep  Aqp 


It  then  follows  from  equations  (5.1),  (5*2),  (5*3) 


Atp 


Atr 


1 

c 


dSp  +  — 
r  c 


(5A) 


This  equation  is  the  analog  to  equation  (4.9) „  From  now  on  the  derivation 
is  very  similar  to  that  in  Chapters  III  and  IV.  The  final  result  for  Case 
III  is 


<$m  =  (  V1  -  Pe  ’)  .  1  -  K  Pm  cos  V 

“ee  <  V  1  -  Pm  '>  1  "  <ne  ^e  cos  *e> 
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The  result  for  Case  IV  is 


wee 


<  V 1  -  Pr  >  "  W  Pm  cos  VJ  “  <ne  Pe  cos  *e>)  /r  ^ 

_  =  - — : — —  ‘  - - - * - Z -  ^-o) 

wrr  <  V  1  “  Pe  >  C1  “  W  Pm  cos  *m) )  f1  “  (nr  Pr  cos  *r) I 


These  equations  replace  the  former  results  (1.9)  and  (1.12).  The  mean 
values  containing  subscript  m  have  already  been  defined  in  Chapters IIII 
and  XV.  Those  containing  subscripts  e  and  r  are  defined  as  follows: 


<  V1  -  Pe  ’>  = 


(t. 


ep 


^eo^ 


^ep  _ 

J  V1  "  Pf(t) '  dt 


eo 


(5.7) 


^rp 

( -Jl  -  )  =  — - 1 -  r  -yji  “  Pr(t)'  dt 

(t.  -  “  '  J 


rP  °ro 


t_)  - 

1X5 '  t 


ro 


(5-8) 


(%  Pm  cos  tn)  = 


(t 


mP 


t  ) 
mo'  -t 


r  mP 

J  %(t)  pm(t)  cos  *m(t)dt  (5.9) 


mo 


tmP 


—  n  air  ^ 

(nm  Pm  COS  tm)  - - — -  /  HjjjCt)  pm(t)  cos  *m(t)dt  (5.10) 

/T  T  \  ^  — 


(’W3  -  Vo^  t 


mo 
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(nr  pr  cos  fr)  =  — — 


(t 


rP 


^rp 


nr(t)  pr(t)  cos  tr(t)dt  (5-11) 


“ro 


Of  course,  it  holds: 


t  =  t 
mo  mo 

tmP  =  tmP. 


The  bars  are  retained,  for  formalistic  fanaticism  only.  However,  they 
prove  quite  useful  in  case  of  t  and  trp,  for  if  we  now  identify  E 
with  R  those  times  become  teo  and  tep  and  these  are  the  times  when 
the  reflected  zero-peak  and  p-peak.  arrive  at  E.  They  are  to  he  distin¬ 
guished  from  teo  and  tgp . 
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CHAPTER  VI 


ACCURACY  CONSIDERATIONS 


So  far  all  formulae  are  exact  (with  the  possible  exception  that  the 
neglect  of  the  boundary  layer  influence  might  result  in  certain  inaccu¬ 
racies).  Paradoxical  as  it  may  sound,  one  starts  to  make  approximations 
as.  soon,  as  one  indulges  in  accuracy  considerations.  We  shall  now  esti¬ 
mate  the  error  in  the  velocity  measurement  caused  by  the  error  in  the 
measurement  of  the  time  of  observation  At. 

We  first  mention  that  the  time  intervals  in  the  definitions  (3.8), 
(4„l4),  (5°7)j  (5 “8),  (5°9)j  (5°10)>  (5-ll)  for  the  respective  mean  val¬ 
ues  are  all  approximately  equal: 

At  =>  tep  -  teo  «  t^  -  tmo  «  trp  -  trQ 

Hence,  all  of  those  mean  values  which  are  basically  mean  values  of 
velocity  components  are  extended  over  approximately  the  same  time  in¬ 
terval  At.  In  other  words,  At  is  the  uncertainty  of  the  time  to 
which  those  velocity  components  refer . 

We  shall  now  investigate  how  the  error  in  the  measurement  of  the 
velocity  components  is  related  to  the  error  in  the  measurement  of  At. 

We  confine  this  consideration  to  Case  II  of  Chapter  I. 

The  basic  assumption  is  that  the  error  dAt  (the  symbol  d  in  this 
chapter  is  used  to  denote  errors)  in  the  measurement  of  the  duration  of 
observation  At  does  not  depend  on  At  itself.  Within  certain  limits 
this  is  in  good  agreement  with  reality. 
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Let  At  denote  the  "true  value"  and  At*  the  "measured  value"  of  the 
duration  of  observation .  We  have  then 

At  =  At*  +  dAt  (6.l) 

Of  course ,  we  assume 


12£*L « i2“L « !  (6.2) 

At  At* 

We  now  apply  this  to  equation  (2.8)  which  holds  for  Case  II  of  Chapter  I. 
This  yields 


K*WI 


2n  g _ 

u>ee  (At*  +  dAt) 


or,  because  of  relation  (6.2) 


1  ^mn^  I 


2  a  p  A  dAt  \ 

uee  At  L  At  J 


for  which  we  may  write 


I  <Pmn>  I 


_  2rt  g 

”  coee  At 


+  ld(Pmn)l 


with 


ld<Pmn)l  - 


2it  a 
we  e  At 


j  dAt  | 

~ At""" 
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If  we  substitute  here  equation  (2.8)  we  obtain 


d(Pmn) 

d<vmn> 

|  dAt| 

^Pmn^ 

^VTnn) 

At 

“eel  fern)  1 
2it  a 


|dAt| 


(6.3) 


This  equation  relates  the  relative  velocity  error  to  the  relative  error 
of  the  time  measurement.  If  we  rearrange  it  a  little  we  also  obtain 


Nvmn)l  At  =  1^1  *  Kvmn>l  (6A) 

In  this  form  the  error  relation  reminds  one  of  Heisenberg's  uncertainty 
relation  which  is,  for  instance. 


AE  At  b  h  (6.5) 

where  AE  and  At  are  the  uncertainties  of  the  energy  E  of  an  ele¬ 
mentary  particle  and  the  time  t  to  which  E  refers . 

The  similarities  between  the  two  equations  (6.4)  and  (6.5)  or  between 
the  corresponding  physical  situations  are  intriguing  even  though  not  un¬ 
limited.  On  the  left-hand  sides  of  both  equations  stand  the  products  of 
two  uncertainties.  In  equation  (6.4)  these  are  the  uncertainties  of  the 
missile  velocity  and  of  the  time  to  which  the  velocity  refers.  On  the 
right-hand  sides  of  both  equations  stand  quantities  which  are  not  unlim¬ 
itedly  at  our  disposal.  However,  it  is  here  where  the  analogy  eventually 
stops.  While  the  natural  constant  h  is  not  at  all  at  our  disposal,  the 
quantity  |dAtj  is  dictated  by  the  technique  of  measurement. 
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The  physical  explanation  for  the  existence  of  the  two  uncertainty 
relations  is  also  rather  analogous.  To  measure  the  energy  of  a  particle 
means  to  measure  its  frequency  according  to  Einstein's  law 


For  measuring  a  frequency  one  needs  a  certain  time  of  observation.  The 
greater  one  selects  this  time  of  observation  the  smaller  becomes  the 
frequency  or  energy  error,  but  the  greater  becomes  the  uncertainty  of 
the  time  to  which  the  energy  refers.  The  situation  is  quite  similar  for 
the  measurement  of  the  missile  velocity  which  now  plays  the  role  of  the 
particle  energy.  Besides,  particle  energy  and  missile  velocity  can  both 
be  translated  into  frequencies . 

This  sheds  a  new  light  on  our  results  (4.15),  (4.16),  (5.5),  and 
(5.6).  The  two  facts  that  ]dAtj  is  dictated  by  the  state  of  the  art 
and  that  we  cannot  tolerate  unlimited  errors  1  d{vmn )  1  mean  that  we 
cannot  measure  the  missile  velocity  as  a  continuous  function  of  time 
but  only  as  a  step  function  consisting  of  a  sequence  of  mean  values. 

This  is  an  inherent  limitation  of  the  whole  method  of  measuring  veloc¬ 
ities  by  means  of  the  Doppler  shift. 

It  is  worthwhile  to  consider  a  numerical  example.  We  shall  vary  a 
few  parameters  but  keep  the  following  parameters  constant  throughout  the 
whole  example: 

The  missile  velocity  vm  shall  be  such  that 


l  <p„)l  - 10-5 

34 


(6.6) 


The  angular  radar  frequency  shall  be 

wee  =  8rt  X  10®  s”1  (6.7) 

This  determines  the  minimum  time  of  observation 


(6.8) 


In  order  to  have  a  completely  determined  situation  we  have  to  make 
two  further  assumptions .  We  do  this  in  two  different  ways .  Once  we  as¬ 
sume  a  fixed  value  for  the  relative  velocity  error 


and  vary  the  time  error  |dAt|,  and  once  we  assume  a  fixed  value  for  the 
duration  of  observation  At  and  vary  again  |&At|  .  In  the  first  case. 
At,  a,  and  P  are  functions  of  the  assumptions,  in  the  latter  case 


d(v  ) 
mn' 


a,  and  P 


is  presented  in  Table  II 


are  functions  of  the  assumptions . 
,  the  second  in  Table  III. 


The  first  case 


35 


TABLE  II 


FIXED  ASSUMPTION 


&{v  ) 

4  mrr 

^vmn) 


=  10 


-6 


|<3At| 

s 

At 

s 

a 

P 

10"6 

10° 

4  x  io4 

4  x  io9 

10"7 

IO"1 

4  x  io3 

4  x  io8 

10“a 

10  “2 

4  x  io2 

4  x  io7 

TABLE  III 
FIXED  ASSUMPTION 

At  =  10”1  s 


3  6 


We  selected  the  three  values  10-6,  10“ 7,  and  10-8  s  for  |dAtj  in 
order  to  satisfy  the  three  psychological  types  of  engineers  -  the  con¬ 
servatives,  the  even-tempered,  and  the  optimists. 
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